Scattering in the nN negative parity channel in lattice QCD 
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We study the coupled irN system (negative parity, isospin |) based on a lattice QCD simulation 
for n/=2 mass degenerate light quarks. Both, standard 3-quarks baryon operators as well as meson- 
baryon (4+l)-quark operators are included. This is an exploratory study for just one lattice size 
and lattice spacing and at a pion mass of 266 MeV. Using the distillation method and variational 
analysis we determine energy levels of the lowest eigenstates. Comparison with the results of simple 
3-quark correlation studies exhibits drastic differences and a new level appears. A clearer picture of 
the negative parity nucleon spectrum emerges. For the parameters of the simulation we may assume 
elastic s-wave scattering and can derive values of the phase shift. 

PACS numbers: 11.15.Ha, 12.38.Gc 



I. INTRODUCTION 

Even if we consider only strong interactions almost all 
hadrons are unstable. Calculations in lattice Quantum 
Chromodynamics (QCD) should therefore take into ac- 
count the resonant nature of these states and the coupled 
decay channels. The bulk of lattice studies rely on cor- 
relation functions for simple qq or qqq-type operators for 
mesons or baryons, respectively. Formally one would ex- 
pect that in the full quantum field theory with dynamical 
quarks these simple meson or baryon operators should 
(via dynamical vacuum loops) couple to meson-meson or 
meson-baryon states. It was found that such intermedi- 
ate channels seem to be coupling too weak to be observed 
(see, e.g., [2-0] for baryon correlation studies). Therefore 
one needs to include explicitly hadron-hadron operators 
in the set of interpolators, as has been demonstrated in 
meson resonance studies |8l4l4| . 

The interplay between resonance levels and hadron- 
hadron states has been discussed in flB - flci ] , where the re- 
sulting energy levels for finite spatial volume were related 
to the continuum scattering phase shift in the elastic re- 
gion. Comparing the energy levels of a non-interacting 
hadron-hadron state with those in the case of interac- 
tions one finds a significant level shift ("avoided level 
crossing") in the energy region of the resonance. The ef- 
fect of such coupled channels depends on the system pa- 
rameters. For small volumes and unphysical large quark 
masses the two-hadron energy levels may lie high above 
the observed resonance levels or - for narrow resonances 
- outside the influence region of the resonance. 

Often it is technically not possible (e.g., due to a small 
volume) to determine more than a few lowest energy lev- 
els below the elastic threshold. In the elastic scatter- 
ing region each energy level corresponds to one value of 
the phase shift and the resonance region then cannot be 
mapped out sufficiently well. One the other hand, each 
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change of volume or other parameters requires a com- 
pletely new simulation sequence (i.e., generating configu- 
rations with dynamical fcrmions, quark propagators etc.) 
Studying interpolators in moving frames [19h23| allows to 
obtain further values on the same configurations. Unfor- 
tunately, for coupled channels with two hadrons of differ- 
ent mass, there can be mixing between different partial 
waves, complicating the situation. Another complication 
is the opening of inelastic channels. 

Starting from continuum models for a scattering pro- 
cess, based on phenomenologically determined parame- 
ters, one can also derive the energy levels on finite volume 
lattices [2^ - |27| . Coupled channel potential models or 
Unitarized Chiral Perturbation Theory motivated mod- 
els in that way allow to compare with lattice results. Al- 
ternative methods to identify resonance parameters have 
been discussed in that context [ID, [28|, Hal ■ 

A particularly interesting case is the negative parity 
nucleon channel. There we have two low lying resonances 
N*(1535) and A^*(1650) which couple to Ntt in s-wavc. 
Above the 10% level there are also further inelastic decays 
A* (1535) -> Nri and AT* (1650) -> Nr],AK. So far lattice 
simulations of this channel, that have determined ground 
state energy levels and further excitations, included only 
3-quark interpolators 043 • I n these studies two low lying 
energy levels have been identified and assigned to the two 
negative parity resonances. However, the lower of the two 
levels showed a tendency to lie below the A*(1535). 

In order to clarify the situation we study here for the 
first time the coupled system of 3-quark nucleon inter- 
polators and pion-nucleon interpolators in the negative 
parity channel. The calculation requires the computation 
of many more correlation graphs than before, including 
the notoriously demanding backtracking quark line con- 
tributions. We therefore use the distillation method [3(| 
for determining the cross correlation matrix for up to 9 
interpolators. We use gauge configurations with rif = 2 
mass degenerate dynamical quarks (of improved Wilson 
type) with a pion mass of 266 MeV. The 16 3 x 32 lat- 
tices have spatial extent of 1.98 fm with Lm, sa 2.68 (for 
details see Table H]) . The energy levels are obtained with 
the variational method I15l [31 -133 [ . 
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Ng x Nt P o[fm] L[fm] Lm^ #configs m^pVleV] 



16 J x 32 7.1 0.1239(13) 1.98 2. 



280 266(3) (3) 



TABLE I: Configurations used for the present study. N$ and 
Nt denote the number of lattice points in spatial and time 
directions, and L — Nga is the size in physical units. 



The energy levels of the eigenstates in case of a finite 
spatial extent L are discrete. They are determined by di- 
agonalizing the correlation matrix of interpolating opera- 
tors. The set of these interpolators has to be large enough 
to allow the representation of the eigenstates. For total 
momentum zero the pion-nucleon operators will have the 
form N(n)n(—n) where n abbreviates the possible quan- 
tized momentum values 2mr/L. For the non- interacting 
situation the corresponding energies are straightforward 
to compute, for the interacting case they are shifted and 
have to be determined numerically from the correlation 
matrix. We need to consider all interpolators that may 
couple to the system in the energy region where one ex- 
pects eigenstates. Obviously the 3-quark interpolators 
and the interpolator iV(0)7r(0) have to be included. In 
our setting already the s-wave operator JV(l)7r(— 1) lies 
high above the lowest energy level. The same holds for a 
possible Nr/2 channel (note, that for only two dynamical 
quarks there is just one rj meson called 772). We find that 
the spectrum shows a clear difference whether the pion- 
nucleon operator is included or not. If the pion-nucleon 
interpolator is included we observe one more level below 
threshold, typical for attractive channels, and the next 
two levels are shifted closer to the expected resonances. 

Following Sect. [TT] where we present the parameters 
and methods used, we discuss the results in Sect. Mil 
The appendix lists the necessary Wick contractions for 
the meson-baryon correlators. 



II. METHODS 

A. Lattice action and configurations 

We use configurations from the study of re-weighting 
techniques [34l |35| generously provided by the authors. 
The gauge configurations were generated for rif = 
2 flavors of mass-degenerate light quarks and a tree 
level improved Wilson-Clover action with gauge links 
smeared using one level of normalized hypercubic smear- 
ing (nHYP smearing). The valence u/d quarks have the 
same mass as the sea u/d quarks. TableUlists the param- 
eters used for the simulation along with the number of 
(approximately independent) gauge configurations used, 
the lattice spacing, volume and the pion mass (for details 
see [H, |Tl|). We note that the small value Lra^ m 2.68 
may lead to finite size effects which we cannot identify in 
this study, since we have just one lattice size available. 



B. Determination of energy levels 

Due to the finitcness of the spatial volume, the en- 
ergy spectrum of the correlation functions is discrete. 
We determine the energy levels of the N and the N it 
system with the variational method [TH, l3ll - [33| . For 
a given quantum channel one measures the Euclidean 
cross-correlation matrix C(t) between several interpola- 
tors, 



djit) = (Oi(t) O 3 -(0)) = l» 



-E„t 



"10,(0)) 



(1) 

where the operators are located on the corresponding Eu- 
clidean time slices. The generalized eigenvalue problem 
C(t)u n (t) = X n (t)C(to)u n (t) disentangles the eigenstates 
\n) with the eigenvalues 

\ n (t,t ) = e- E ^-^(l + o(c- AE ^-^)) , (2) 

where AE may be as small as the distance to the next 
nearby energy level. From the exponential decay one 
determines the energy values of the eigenstates by expo- 
nential fits. The stability of the eigenvectors with regard 
to t and the so-called effective energies 



E n {t) = log 



\ n (t+l) 



(3) 



indicate the suitable fit range by exhibiting plateau-like 
behavior. The set of interpolators should be large enough 
to allow the system to reproduce the physical eigenstates. 
Neglecting important interpolators may obscure the re- 
sult. On the other hand, in the calculations it is not pos- 
sible to have a complete set of interpolators and one is 
limited to a reasonable subset. Also, the statistical qual- 
ity of C(t) is an issue. The reliability of the obtained 
energy levels decreases for higher \n), with the ground 
state being the most reliable one. 

The energy values are extracted using correlated fits of 
X n (t) to one and two exponentials. A possible source of 
systematic error is the choice of the fit range in t. From 
the effective energy plots (cf., Figs. [2] and [3]), the range of 
stability of the eigenvectors and the \ 2 dependence of the 
fits we estimate suitable fit ranges. The two exponential 
fits start at smaller t and we verify that the extracted 
levels agree with results obtained from one-exponential 
fits starting at larger t. 



C. Interpolators 



The Nir system can be projected to isospin i and |, 
experimentally accessible by ir^p scattering. Here we 
study only the isospin i sector. 

For the charged nucleon interpolator we use the opcr- 
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ator (on a given time slice) 
(JV«) M (p = 0) = 

5> 6c (p±r« Uo (2)) (^(f)r^d c (f)) , (4) 

5? 

and for the neutral one with the quarks d, u, d. (Ti, I^) 
can assume the three values (1,6*75), (75, C) and 
(il,C7 f 7 5 ) for i = 1,2,3. C denotes the charge con- 
jugation matrix, 74 the Dirac matrix in time direction, 
and P± = ^(1 ± 7 t ) the parity projector. We sum over 
all points of the time slice in order to project to zero 
momentum. Summation over the color indices a, 6, c and 
the (not shown) Dirac indices is implied. 

In the distillation approach (see Subsect. Ill Dl below) 
the sources are smeared combining N v eigenvectors. For 
the nucleon 3-quark interpolators we choose N v = 32 
and N v = 64 and thus with the three different Dirac 
structures have six operators. 

The pion interpolators read 

tt + (p = 0) =^2d a (x)j 5 u a (x) , 

X 

ir°{p = fy=y~\-7= {u a (x)i 5 u a (x) - d a (x)j 5 d a (x)) , 

(5) 

where summation over the color index a is implied. 

We consider the Nir system in the rest frame. The 
leading s-wave contribution then comes from the inter- 
polator with both particles at rest, 

Nir(p = Q) = 75 jV + (p = 0)7r(p = 0) , (6) 

where N+ denotes the positive parity nucleon and the 
factor 75 ensures negative parity for the interpolator. In 
the distillation approach we choose for the Nir channel 
N v = 32 and thus with the three different nucleon inter- 
polators have three operators. 

We project to isospin i by choosing the combination 

OnM = \' I * = \)= + ^ n7T+ > ( ? ) 

with p and n denoting the charged and the neutral nu- 
cleon according to (J4]). 

The negative parity TV* channel becomes quickly in- 
clastic (see, e.g., |3a - t39| ). According to the Particle Data 
Group [4(| the main decay channel is Nit (35-55% for 
iV*(1535), 50-90 % for JV*(1650)). The second largest 
decay rate is to Ntj ( 42± 10% for JV*(1535), 5-15 % for 
JV*(1650)). Most of the rest of 10-20 % is Ntttt and, for 
iV*(1650), also AK. For a lattice calculation at physi- 
cal quark masses one would need to include the inelastic 
channels. This is beyond present days capacities. 

In our case (n/ = 2) there is just one pseudoscalar 
meson 77 called 772 , with a mass larger than 800 MeV [U, 
|42| . With our parameters (see Sect. IIIip these inelastic 



channels would thus have thresholds above 1900 MeV. 
The lowest state with total momentum zero but non- 
zero relative momentum iV(l)7r(— 1) (momentum units 
2-7r/(16a)) has a (non-interacting) energy of 1920 MeV, as 
well. These energy values are above the observed lowest 
three levels. We cannot exclude that in particular the 
highest of these may be influenced by the N(1)tt(—1) 
state. 



D. Distillation method and correlation function 

We compute the correlation matrix entries with help 
of the distillation method j3(j. This method has been 
successfully applied in several studies, including baryon 
correlation functions @, II @, [H [H El 13 . It also al- 
lows for a reliable evaluation of the partially disconnected 
diagrams. On a given time slice one introduces separa- 
ble (i.e., expressed by a sum of products separating the 
dependence on x and x') quark smearing sources in the 
form 

Qc,a {x ) -> ^ S c d(x, x')qd, a {x') 
x' 

N v 

^EE^^w*'), (8) 

x' i 

where c, d and a denote color and Dirac indices and sum- 
mation over the color indices is implied. A suitable choice 
for the v % is the eigenvectors of the spatial lattice Lapla- 
cian [30( | . Summing over all eigenvectors reproduces the 
delta function, the spectral representation of unity. In 
actual calculation one truncates the sum and uses the 
lowest eigenmodes or subsets. The value of N v depends 
on the lattice size and N v between 32 and 96 was found 
suitable for our situation Q. 

The advantage of the distillation approach lies in its 
versatility. Instead of quark propagators G CK d v {x,Xo) 
from one source located in xq to other points on the lat- 
tice one now computes propagators between eigenmode 
sources, so-called perambulators 

7~fj.v(.j: ^snki ^ 7 ^-sre) — 

E v li{^^^\dGd^cu{x,y^v l c {y,t^ c ). (9) 

The interpolator structure decouples from the calculation 
of the perambulators completely. E.g., meson correlators 
assume the form 

C(7j s „k,4rc) = (M(t 8 nk)M t (t src )) 

$liv (tTj , k'i ^src)^a ^sre ; ^'j ^snk) 

<l>ap(h j',tsnk)rpn(J, t sri k;ri, t src ) , (10) 

where M denotes a meson interpolator like, e.g., the pion 
of Eq. ([5]) and summation over the source index (i, j, k, n) 
pairs and the Dirac index (a, f3, fx,v) pairs is implied. 
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Due to 75-hermiticity of the Dirac operator the peram- 
bulator for sink to source can be expressed by that from 
source to sink, 

(11) 

The meson interpolator type is specified in 



D a p <i>(i,j;t) 



(12) 



The factors D and F represent the Dirac structure and 
momentum projection or derivative terms related to the 
quantum numbers of the meson. Only <p has to be re- 
computed for each meson interpolator whereas the per- 
ambulator remains the same. 

For 3-quark interpolators like the baryons one obtains 
contributions to the correlation function of the form 

C^(^nk, t SIC ) = (-/V A1 (t sn k)A r l/ (t src )) 
= 'Pp.aflj (i 1 Ji k\ t sn k) 

T aa ' (z, ^snki ^ 5 ^sre) 

Tp/3<(j,t sn k;j ,t SIC ) 
^Lp^(^j\k';t src ) . (13) 



For an interpolator N (without derivatives) 
the form 



assumes 



<l>vap~/{i,j, k] t) = D vafil V e a bcV z a {x)vl(x)Vc(x)F(x) 



DuapjHhj: k; t) 



(14) 



Again, D carries the Dirac structure and F the possible 
total momentum projection factors. 

We also project the correlation functions to definite 
parity with the projection operators P± = j(l ± 7j). In 
App. [A]we list the necessary contraction terms expressed 
in terms of the perambulators. 



E. Energy levels: interpretation 



We study the Nir system in the i channel in s-wave 
in the rest frame. From the energy value 



E = sfs = ^(pn+Pk) 2 = vV 2 + ml + yp* 2 + m% 

(15) 

we extract the momentum p* = |p*| with 

» 2 [s - (m N + m„) 2 ][s - (m N - m n ) 2 } 
P = ^ , (16) 

and the dimensionless product of the momentum and the 
spatial lattice size 




FIG. 1: Lowest energy levels vs. the spatial lattice size (both 
in units of the pion mass). Upper plot: physical pion, nu- 
cleon and TV* (1535) masses, comparing the non-interacting 
levels (dotted) with the levels distorted due to interaction 
(full lines); the broken horizontal line indicates the thresh- 
old, the horizontal thick line the N* (1535)-mass; the N* is 
parametrized as an elastic resonance with a decay width of 
150 MeV. Lower plot: Unphysical values for = 266 MeV, 
mjv = 1068 MeV; for TV* the mass is chosen as 1670 MeV 
without changing the coupling strength. In both cases the 
lowest possible state is 7V(0)7r(0), which coincides with the 
threshold in the non-interacting case. For attractive interac- 
tion the level moves slightly below the threshold to negative 
2 



q=p 



L 
2^ 



(17) 



For a system of non-interacting pions and nucleons 
the energy levels for given lattice size can be straight- 
forwardly computed (dotted lines in Fig. [1]). For the in- 
teracting case, with localized interaction region and in 
the elastic domain, Luscher fl5l - [l8j has given a relation 
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between energy levels and phase shift, 
tan S(q) = 



ir 3 / 2 q 



2bo(i; q 2 ) 



(18) 



where q is given in (jTTJ) , and the generalized zeta function 
Z\ m is given in (HI- 

Assuming a phase shift parameterization, one can nu- 
merically invert that relation and obtain the modified 
energy levels, which exhibit the phenomenon of avoided 
level crossing by level "transmutation" . 

The s-wave amplitude may be written 



T = c ld sin 8 = 



1 



cot S — i 



We also define for convenience 



p* 2Z 00 (l; g 2 
po(s) = —= cot 5 = — / — 



(19) 



(20) 



with the effective range parameterization near threshold 

1 



1po(a) = — + 0(p*'O , 
ao 



(21) 



and scattering length ao . If the first resonance is of Brcit- 
Wigner shape, then po can be approximated linearly, 



Po{s) = ~(sr - s) 

7 



(22) 



Here sr denotes the resonance position and 7 is related 
to the width 



P*(sr) 



sr 



(23) 



or the coupling constant 7 = g 2 /6ir. 

The iV(g ) (s-wave) scattering amplitude is shown in 
the data analysis of |36[ and has an intricate behavior, 
becoming quickly inelastic. In a simplification of that 
case let us study the situation with just one elastic reso- 
nance. In that case the phase shift and elastic amplitude 
can be modeled where we have used a resonance mass of 
1535 MeV and a width of 150 McV. The resulting energy 
levels demonstrating the expected avoided level crossing 
are also shown in Fig. Q] 

In the lower part of Fig. [I] we show the situation where 
the pion mass has the larger value 266 MeV. (Note that 
it is also used as unit mass in that plot.) The values 
of the stable nucleon has been set to 1068 MeV and the 
resonance position to 1670 MeV, all values close to the 
results of our calculation to be discussed in Sect. IIIII 
The coupling strength 7 at the resonance position is un- 
changed. 

In this setting the picture changes drastically. For our 
lattice size we have L m 2.68 (in units of m ff ) and the 
energy level iV(l)7r(— 1) lies clearly above the resonance. 
The lowest possible state is A(0)7r(0), which coincides 



with the threshold in the non-interacting case. For at- 
tractive interaction the level moves slightly below the 
threshold to negative q 2 , which is a finite volume arti- 
fact. 

Choosing interpolators with non-zero total momentum 
( "moving frame" ) allows in principle to obtain further en- 
ergy levels and thus additional values of the phase shift. 
For the case of two particle of equal mass this was dis- 
cussed in 0, [20| and has been used in various studies of 
the 7T7r-system. The situation for pairs of hadrons with 
different masses is more complicated [2lT - 423| since there 
even and odd partial waves may mix. In this study we 
rely on the case of zero momentum. Smaller quark masses 
will require the consideration of further Nt: operators 
and other interpolators. 

Liischer's relation holds in the elastic region. Most 
often inelasticity sets in early due to coupled channels. 
An alternative approach is the inverse procedure, starting 
with a (unitarized) coupled channel parameterization of 
the scattering matrix in continuum and then determining 
the expected discrete energy levels on finite volumes, see, 
e.g., I24I27I ]. The lattice results for the energy levels can 
then be interpreted along these lines. 



III. RESULTS 

A. Pion and nucleon, non-interacting 

The masses of the free pion and the ground state nu- 
cleon A(i + ) have to be estimated with the highest possi- 
ble precision in order to perform the subsequent analysis. 
For the gauge configurations used here the pion was stud- 
ied carefully in [1, EH where the value am^ = 0.1673(2) 
was obtained and we use this value here as well. 



For the positive parity nucleon N(h 
correlation matrix for the six operators 



we study the 



0\,Ol,0\ = vj 1} , N { +\ A| 3) with N v = 32 , 



(3) 



with N„ 



64 



(24) 



The correlation matrix is analyzed with the variational 
method as discussed above. The ground state shows a 
stable plateau behavior in the effective energy plot Fig. [2] 
The first excitation is considerably higher than the ex- 
pected Roper resonance. This observation is share d by 
other recent studies (see, e.g., [3,IH) but disputed [2ll45|. 
The reason for the high value may lie in the incomplete- 
ness of the interpolator basis, i.e., possibly missing im- 
portant 5-quark interpolators. To solve this puzzle is not 
the object of our study. Our value of the ground state nu- 
cleon (fit range 6-12) is am^ — 0.672(4) (corresponding 
to m N = 1068(6) McV). 
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FIG. 2: The effective energy values for the JV(i + ) channel 
(with 3-quark interpolators). 

B. Interacting iV7r system 

We compute the full correlation matrix for the follow- 
ing operators: 



Of, 02 ,0s 

or, , o 6 - 



N W tN M iN W 
iV (1) , iV (2) , iV (3) 



with iV,, = 32 
with N v = 64 
Of, 8 ", 0f = 0<£, 0<£, with iV„ = 32 , 



(25) 

with the definition from (J3J and ((BJ). 

Let us first consider results for the subset of 3-quark 
interpolators Of-Of . It turns out that inclusion of 

(3) 

the type N_ ' does not improve the quality of the di- 
agonalization results. We therefore use only the subset 
(01, 0f , Of , Of)- We reproduce the usual (see, e.g. 
0,13) pattern of energy levels (see left hand plot of 
Fig. [3]), which have been assigned to the two N* reso- 
nances. However, as has been observed in [4o| . towards 
smaller pion masses the lower level moves close to the ex- 
pected threshold and thus lies unexpectedly low if com- 
pared to the JV*(1535). The situation is shown in Fig. 0] 
(middle). The energy levels have the values 1.359(43) 
GeV (exponential fit, fit range 6-10) and 1.709(29) GeV 
(fit range 4-9). 

This picture changes significantly, when one includes 
the iV7r-interpolators in the correlation matrix. The right 
hand plot of Fig. [3] shows the effective energy levels 
when using operators Of , Of , Of , Of , Of , Of , Of in 
the analysis. The exponential fits to the corresponding 
eigenvalues and the resulting energy levels are listed in 
Table HI 

Figurc@] (right) demonstrates the difference to the pre- 
vious case with only 3-quark interpolators. The lowest 
level now lies slightly below threshold, a feature typical 
for attractive s-wave (Til [l2[ and a finite volume arti- 
fact. This agrees with the behavior discussed in Subsect. 
Ill El The next-higher two levels are now close to values ly- 
ing approximately 130 MeV above the physical resonance 
positions of iV*(1535) and iV*(1650), similar to the situ- 
ation for the nuclcon. Comparison with Fig. [TJ where a 
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FIG. 3: Left: effective energy values for the case without Nn 
contribution, right: including Nn interpolators. The horizon- 
tal broken line indicates the threshold value mjy + "V- 




N_ Ntt 



FIG. 4: Comparison of the energy levels. Left: physical mass 
values (experiment). Middle: result when using only 3-quark 
interpolators. Right: result when pion-nucleon interpolators 
are included. The dashed lines indicate the scattering thresh- 
olds. 



single clastic resonance parameterization has been used, 
shows excellent agreement for the lowest two energy lev- 
els. 

The eigenvectors are fingerprints of the states and one 
should have a stable composition across the fit range in 
order to be sure to identify the same eigenstate. Fig. [5] 
shows the eigenvector components of the three lowest 
eigenstates. The eigenvectors have unit norm. The ab- 
solute normalization of the 5-quark operators compared 
to the 3-quark ones is unclear. However, one finds that 
the contribution to the ground state is significantly 
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level 
n 


to 


fit 
range 


a E n = a sj~s 


E = </s 
[GeV] 


x 2 

d.o.f. 


ap* 


Po 


S 

[degrees] 


1 


1 


6-12 


0.800(5) 


1.272(8) 


6.12/5 


0.0985(57) i 


0.149(48) 


68(59)i 


2 


1 


4-8 


1.045(19) 


1.662(30) 


2.46/3 


0.2726(155) 


0.007(42) 


89(9) 


3 


1 


4-8 


1.127(18) 


1.792(29) 


0.67/3 


0.3362(134) 


0.279(108) 


47(10) 



TABLE II: Final results for the lowest three energy levels of the coupled Ntt system with the interpolators 
Oi ,C>2 ,0~[ ,0$ ,Oi ,Og ,Oq . The energy levels are determined by correlated one-exponential fits to the eigenvalues X„(t) in 
the given fit range. We verified that two-exponential fits starting at smaller t agree with results obtained from one-exponential 
fits. The errors are determined by the single-elimination jackknife method. For the values given in GeV we use the lattice 
spacing a — 0.1239 fm (Sommer parameter ro = 0.48fm). 



larger than to the higher levels. Interpolators of type 
contribute importantly to the lowest eigenstate and 
dominate the 3rd state, whereas the interpolators of type 

(2) 

N_ ' are more important for the 2nd state. 

In contrast, the effective energy levels of the pure 3- 
quark correlations system show more fluctuation. Com- 
paring with the full Ntt system results one gets the im- 
pression that the two lowest states of the 3-quark system 
interpolate between the three lowest states of the com- 
plete system. 

The lowest energy level of the two particle system lies 
below threshold and the corresponding value of po may 
be related to the scattering length. Table|TT]gives also the 
values of po from (|2TJ| due to the Luscher analysis and 
the resulting values of the phase shift, assuming elasticity. 
The second energy level lies close to the point where the 
phase shift crosses ir/2 (this value is included within the 
error bars). This closeness is pure chance: for slightly 
larger lattices this would not have been the case (cf., 
Fig. [1}. As discussed, the kinematical situation (pion 
mass and lattice size) allows the assumption to be in the 
elastic domain and thus one is tempted to assume validity 
of (|22l) . The zero of the line connecting the values of po at 
the two lowest energy levels give the resonance position 
cl 2 sr = 1.114(135) corresponding to a resonance mass 
m R = 1.678(99) GeV. This is approximately 140 McV 
above the physical value, but not surprising due to the 
unphysical pion and nucleon masses, in fact, a similar 
shift as for the nucleon. Also note, that the Ntt system in 
Nature is already inelastic and the linearity assumption 
not justified in that case. 

The third eigenstate has a phase shift of 47° (~ 227°, 
since the arctan is defined modulo 180°), indicating a 
resonance lying closely above that energy value of 1.79 
GeV - again assuming elastic scattering. 

Due to the closeness of the threshold to the resonance 
in our setting, as compared to Nature, we cannot ex- 
pect physical values for scattering length or decay width. 
With (|2"Tj) we can estimate the scattering length from 
the point close below threshold s t /, r . We find a value 
a ~ 5. 3 (±1.4) GeV -1 roughly four times larger than, 
e.g., the lea ding order Chiral Perturbation Theory value 
m ff /(47rF2) [43E3. 



IV. SUMMARY 

We studied Ntt scattering in the negative parity, 
isospin i sector in an ab initio lattice QCD calculation. 
The simulation parameters are: two dynamical, mass de- 
generate quarks, a pion mass of 266 McV, a spatial lattice 
size of 1.98 fm, a volume 16 3 x 32 in lattice units. We 
use 3-quark and meson-baryon (5 quark) interpolators 
and analyze the 9x9 correlation matrix with help of the 
variational method. 

We find a significant difference to the results of sim- 
ple 3-quark correlation analyses. The overall behavior is 
resembling that found in meson-meson scattering lattice 
studies for s-wave channels [ill, 113 • Due to the unphys- 
ical values of the pion mass, the resonance position is 
higher than the experimentally established values. 

The main result of our study is that taking into account 
meson-baryon interpolators indeed changes the obtained 
energy spectrum significantly. This is a first step into 
that direction. Obviously this is an exploratory study 
and systematic uncertainties stemming from the volume 
size, the lattice spacing and the pion mass are not (yet) 
under control. More work (moving frames, different vol- 
umes, further coupled channels) will fill the gap between 
clastic and inelastic threshold and allow the comparison 
with experiment and continuum models. 
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reads 



FIG. 5: The (normalized) eigenvector components for the low- 
est three eigenstates observed; the f-range used for the expo- 
nential fit to the eigenvalues is indicated by a broken line. In 
the legends the operator numbers according to (|25|) are given. 



Appendix A: Wick contractions 

Notation for the perambulators used in this sec- 
tion: T(t,t',a,a',a,a') denotes the perambulator 
T aa '(a,t;a',t') from ©, i.e., from source at t' (source 
vector a', Dirac index a') to the sink at t (source vector 
a, Dirac index a). 

Each source/sink nucleon contributes a factor of the 
form ^Ar(a, b, c), which is constructed from the Laplacian 
eigenvectors. For a given time slice we have 



<i>T{a,b,c) 



k vl{x)v 3 b {x)v k c {x) 



(Al) 



where e denotes the Levi-Civita symbol, v are the Lapla- 
cian eigenvectors, and the sum runs over all sites of the 
time slice and over the color indices i,j,k. The corre- 
sponding factor for the pion 4> n (a, b) on a given time slice 



^ k (a,b) = ^S ij vi*(g)vl(x). (A2) 



By permuting and renaming the Dirac indices a, /3, 7, . . . 
and the eigenvector indices a, b, c, . . . we group the differ- 
ent contraction such that they have a common prefactor. 
There also the gamma matrices of the nucleon and pion 
and the parity projection operators are located. 



1. N -> N 



This entry has the form 



r^ M ^rlrf 7 r^^r (a, b, c)^V, 6', c') £ A t 



where summation over index pairs is implied. 



! = 1 



(A3) 



Source Al 



So u rc e A2 



FIG. 6: Terms Ai and A 2 contributing to N N. 



T(t,t',b,a',l3,a') 
A2 = —t (i, t , a, a', a, a) t (t, t , b, b , ft, (3') 

r(<,t',c, C ', 7 ,7') (A4) 



2. N -> TVtt 



This matrix element has 4 terms contributing: 



72 



p^-f p± pA -p_B p-Bf px 



isnk i 
N 



\a,b,c)ct>T k (e,g)^ c (a',b',c') , (A5) 

i=l 

where summation over index pairs is implied. 
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FIG. 7: Terms Bi - B4. contributing to N -t Ntt. 



B\ = 3 t (t, t, b, e, j3, e) r (t, i', a, c', a, 7') 

T(t,t?,g,a',6,a') T{t,t',c,b',j,0') 
B 2 = -3t (t, t', b, b\ (3, /?') t (t, t', c, c', 7, 7') 

t (t, t, a, e, a, e) r (t, i', g, a', (5, a') 
S 3 = -3r (i, t, c, e, 7, e) r (i, t', b, a', /3, a') 

r(i,t',a,c',a,y)T({,f',j,i',«S,i3') 
Bi = 3 r (t, c, c', 7, 7') t (t, t, a, e, a, e) 

T(t,t',b,a',/3,a') T(t,t',g,b',5,p') (A6) 




FIG. 8: Terms & - C4 contributing to Ntt -t N. 



d = 3 r (t, t', b, b', P, /?') r (t, t', c, c', 7, 7') 

r i', e', a', e', a') r (i, i', a, 5', a, (5') 
C 2 = -3 r (t, i', c, c', 7, 7') r (t, t', a, b', a, (3') 

T(t',t',e',a',e',a') r(t,t',b,g',l3,8') 
C3 = — 3 r (t, i', c, a', 7, a') r (i, t', a, g' ', a, S') 

T(t,t',b,c'^ n ')T(t',t',e , ,b , ,e',^) 
C4 = 3 t (i, i , a, 6', a, /3 ) r (t, t', c, a', 7, a') 

T{tJAsf,W) T(t',t',e',c',e , ,j') (A8) 



3. iW — >■ TV 



4. iW -> JW 



This matrix element has 4 terms contributing: 



Here 19 terms contribute: 



v/2 



r Af p ± r A r B r Bt pirf 
- 1 - ~' ju^ fiv L va L fii L y'0' e'S' 



#?*(M,<#£ c (a' >' ,c')0r(3>e) K» ( A7 ) 



_p^t p± pA pB p-Bf p-n- pwf 
1 /37 1 T'/S' 1 <5 £l 



19 



(a, 6, c)&» fc (e, (a', c'WJ c (g\ e') £ A 



(A9) 



where summation over index pairs is implied. 



where summation over index pairs is implied. 
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FIG. 9: Terms Di - D 5 contributing to Nn -t Nn. 



D 1 =3T(t,t',b,b',/3,/3') T(t,t',c,c',w') 

t (*', t, e', e, e', e) r (t, t' , a, g', a, 5') 

t (t, t' , g, a' , 5, a') 
D 2 = -3t (t 1 , t, e', e, e', e) r (t, t', a, &', a, /?') 

T (t, t', g, a', S, a') r (t, t' , b, c' , /3, 7') 

T{t,t',c,g , ,j,S') 
D3 = — 3 t (t 1 , t, e', e, e' , e) r (t, i , c, a', 7, a') 

D 4 = 9 t (t', i, e', e, e', e) r (i, i', a, c', a, 7') 

rfoi'.c, a',7 ia ')T(M',M',/3^') 
T(t,t',g,b',5,p') 
D5 = — 6 r (£, t', a, a', a, a') t (i', t, e! , e, e', e) 
T(t,t' 1 b 1 c , ,p, 1 ')T(t,t',g,b , ,6,(3') 
T(t,t',c,g', 7 ,S') (A10) 




FIG. 10: Terms D 6 - D 10 contributing to Nn iVvr 



£> 6 = -6 r (f, i, e', e, e', e) r (t, i', 3 , </, S, 5') 
t (t, t', b, a', /3, a') r (i, t' , a, c', a, 7') 

t(M',c,&',7,/3') 
L> 7 = 6 t (i, t', a, a', a, a') r (t, t', b, 6', /3, (3') 

r (i, i', c, c', 7, 7') t (t\ t, e', e, e, e) 

T(t,t',g,g',S,6') 
Dg, = — 9 t (£, t, c, e, 7, e) r (t, t', a, c', a, 7') 

t (t\ t', e', a', e', a') r (i, t', b, g' , (3, S') 

r(t,t',g,b',d,j3') 
Dg = 9 r (t, t', c, c', 7, 7') t (t, f , a, e, a, e) 

T(^^^^e^a^ e ^a')T(M^M^/M') 

T{t,t',g,b',5,p') 
D 10 = 9t (t, t', g, g', S, 5') t (t, t, b, e, 0, e) 

t (t, t', a, c', a, 7') t (t', t', e', a', e', a') 
r(^t',c,6',7,/?') (All) 
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D u = -9t (t, t', b, b', p, 13') t (t, t', c, c', 7, 7 ') 

t (t, t',g, g , 8, 8') r (t, t, a, e, a, e) 

t [t ,t ,e ,a ,e ,a ) 
Du = -3 t (t, t, c, e, 7, e) r (t, t', a, g , a, 8') 

t (t, t',g, a', 8, a') r (t, t', b, c , /3, 7') 

r(t , ,t , ,e , ,b , ,e',p') 
D13 = 3r (i, t, 6, e, /3, e) r (t, i', a, c', a, 7') 

T({,/,D,tt',{, a ')T(t't'e>'^') 

D14 = 3 r (t, i, c, e, 7, e) r (t, t', b, a', /3, a) 

T{t,t',a,g',a,5')T(t',t',e',b',e',p') 
T{t,e,g,<S,6,T') 

D15 = 6 r (t, i', a, a', a, a') t (t, i, c, e, 7, e) 

T(t',t',e f ,b',e',p')T (t,t',b, g',p,8') 
T(t,t',g,c',8,y) (A12) 



Die = -3 t (t, t', g, g', 8, 8') t (t, t, c, e, 7, e) 
t (t, t' , &, a', (3, a 1 ) t (t, t' , a, c', a, 7') 



T(t',t',e',b',e',p') 



D17 = —6 r (i, t' , a, a', a, a') r (i, i', g, g' ', i5, 5') 
T(t,t,c,e,-f,e)T(t,t / ,b,c',/3,^') 
T{t',t',e',b',e',P') 

Dis = 9 r (i, £', g, 5', 5, <5') r (t, i, a, e, a, e) 

r&t'.c, a',^,a')T(t,t',b,c',p,y) 
r(t',t',e',b',e',l3') 

Dig = — 9 r (t, i, a, e, a, e) r (i, i', c, a', 7, a') 
T(t,t',b,g',l3,6')T(t,t',g,b',6,P') 
T(i',t',e',c',e',7') (A13) 
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